EPJ manuscript No. 

(will be inserted by the editor) 



NO 
O 
O 
(N 

-t— > 
O 

O 

ON 



Heavy Quarkonium Physics from Effective Field Theories 

Antonio Vairo a 

Dipartimento di Fisica dell'Universita di Milano and INFN, via Celoria 16, 20133 Milano, Italy 
Received: 19 October 2006 

Abstract. I review recent progress in heavy quarkonium physics from an effective field theory perspective. 
In this unifying framework, I discuss advances in perturbative calculations for low-lying quarkonium ob- 
servables and in lattice calculations for high-lying ones, and progress and lasting puzzles in quarkonium 
production. 
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1 Introduction 

Heavy quarkonium is most frequently and successfully stu- 
died as a non-relativistic bound state of QCD [Q. Non- 
relativistic bound states made of heavy quarks are char- 
acterized by a hierarchy of energy scales: m, mv, mv 2 , 
m being the heavy quark mass and v the heavy quark ve- 
locity in the centre-of-mass frame. A way to disentangle 
rigorously these scales is by substituting QCD, scale by 
scale, with simpler but equivalent Effective Field Theories 
(EFTs). Modes of energy and momentum of order m may 
be integrated out from QCD in a perturbative manner 
(m ^> Aqqd , by definition of heavy quark) leading to an 
EFT known as non relativistic QCD (NRQCD) [20; in- 
tegrating out gluons of energy or momentum of order mv 
leads to an EFT known as potential NRQCD (pNRQCD) 
@|, see Fig.^ pNRQCD is close to a Schrodinger-like de- 
scription of the bound system and, therefore, as simple. 
The bulk of the interaction is carried by potential-like 
terms, but non-potential interactions, associated with the 
propagation of low-energy dynamical degrees of freedom, 
are generally present as well. For a review on the subject 
we refer to EFTs for quarkonium were discussed at 
the conference by J. Soto. 

QCD is characterized by an intrinsic scale ^Iqcd- Let 
us consider states below threshold. In the perspective of 
non-relativistic EFTs, it may be important to distinguish 
between low-lying quarkonium resonances, for which we 
assume mv 2 > /Iqcd, and high- lying states, for which we 
assume mv ~ Aqqd ■ In the first case, the system is weakly 
coupled and the potential is perturbative, in the second 
case, the system is strongly coupled and the potential 
must be determined non-perturbatively, for instance, on 
the lattice. This difference is not appreciated in potential 
models that typically describe the whole spectrum with 
the same interaction. Besides potential terms, which en- 
code the dynamics at the soft scale mv, EFTs contain 
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Fig. 1. EFTs for quarkonium. 



also terms that describe the dynamics of lower energy de- 
grees of freedom. These terms are typically absent in po- 
tential models. For what concerns systems close or above 
the open flavor threshold, a complete and satisfactory un- 
derstanding of their dynamics has not yet been achieved. 
Therefore, the study of these systems, despite their phe- 
nomenological interest, is on a less secure ground than the 
study of states below threshold. Indeed, it largely relies on 
phcnomenological models. 

In the following, I will outline some recent (and very 
recent) progress in our theoretical understanding of heavy 
quarkonium physics. The perspective is quite personal, 
not only because the framework, which I will privilege, 
is that one of QCD EFTs, but also because I can only 
present a selection of results. It is important, when dis- 
cussing about progress, to keep in mind the distinction 
made above about low-lying, high-lying and threshold sta- 
tes. In the first case, we may rely on perturbation theory: 
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the challenge is in performing higher-order calculations 
and the goal is precision physics. In the second case, we 
have to rely on non-perturbative methods: the challenge is 
in providing a consistent framework where to perform lat- 
tice calculations and progress is measured by the progress 
in lattice computations. In the third case, one of the ma- 
jor challenges is to interpret the new charmonium states 
discovered at the B-factories in the last few years (for 
an updated list, see http://www.qwg.to.infn.it/). The dis- 
tinction may help to better appreciate the quality of the 
progress made. 



2 Low-lying QQ 

Low-lying QQ states are assumed to realize the hierarchy: 
m 3> mv 3> mv 2 > vIqcd, where mv is the typical scale 
of the inverse distance between the heavy quark and an- 
tiquark and mv 2 the typical scale of the binding energy. 
At a scale /z such that mv ^> /.i ^S> mv 2 the effective de- 
grees of freedom are QQ states (in color singlet and octet 
configurations), low-energy gluons and light quarks. 

The lowest-lying quarkonium states are r]b (not yet de- 
tected), Y(1S), T] c and J/tp. The T(1S) and J/tjj masses 
may be used to determine the bottom and charm quark 
masses. These determinations are competitive with other 
ones (for the b mass see e.g. |Hj)- We report some recent 
determinations in Tab. 1. 



reference 


order 


m b (rn b ) (GeV) 





NNNLO* 


4.210 ±0.090 ±0.025 





NNLO +charm 


4.190 ±0.020 ±0.025 


ma 


NNLO 


4.24 ±0.10 


m 


NNNLO* 


4.346 ± 0.070 


\m 


NNNLO* 


4.20 ± 0.04 


pi 


NNNLO* 


4.241 ± 0.070 


13 


NNLL* 


4.19 ±0.06 


reference 


order 


m c (m c ) (GeV) 


D31 


NNLO 


1.24 ±0.020 


ma 


NNLO 


1.19 ± 0.11 



Table 1. Different recent determinations of m^ms) and 
rn c (rn c ) in the MS scheme from the bottomonium and the 
charmonium systems. The displayed results either use direct 
determinations or non-relativistic sum rules. Here and in the 
text, the * indicates that the theoretical input is only partially 
complete at that order. 



Once the heavy quark masses are known, one may 
use them to extract other quarkonium ground-state ob- 
servables. At NNLO the B c mass was calculated in ^3] 
{M Bc = 6326 ± 29 MeV), [H] (M Bc = 6324 ± 22 MeV) 
and (M Bc = 6307 ± 17 MeV).' These values agree 
well with the unquenched lattice determination of |16| 
(M Bc = 6304±12tJ 8 MeV), which shows that the B c mass 
is not very sensitive to non-perturbative effects. This is 
confirmed by a recent measurement of the B c in the chan- 
nel B c — > J/t/jir by the CDF collaboration at the Teva- 
tron; they obtain with 360 pb" 1 of data M Bc = 6285.7 ± 



5.3 ± 1.2 MeV \Tj\, while the latest available figure based 
onl.lfb" 1 of data is M Br . = 6276.5 ± 4.0 ±2.7 MeV (see 
http://www-cdf.fnal.gov/physics/new/bottom/060525.ble 
ssed-bc-mass/). 

The bottomonium (and charmonium) ground-state hy- 
perfine splitting has been calculated at NLL in ^SJ. Com- 
bining it with the measured T(1S) mass, this determina- 
tion provides a quite precise prediction for the rjb mass: 
M ni = 9421 ± 10lg MeV, where the first error is an esti- 
mate of the theoretical uncertainty and the second one 
reflects the uncertainty in a s . Note that the discovery 
of the rjb may provide a very competitive source of a s 
at the bottom mass scale with a projected error at the 
Mz scale of about 0.003. Similarly, in ^U], the hyper- 
fine splitting of the B c was calculated at NLL accuracy: 
M B * - M Be = 65 ± 24±ll MeV. 

The ratio of electromagnetic decay widths was calcu- 
lated for the ground state of charmonium and bottomo- 
nium at NNLL order in [201 - I n particular, they report: 
r(r) b -> 77)/r(r(l£0 -> e+e-) = 0.502 ± 0.068 ± 0.014, 
which is a very stable result with respect to scale variation. 
A partial NNLL* order analysis of the absolute width of 
T(1S) — > e + e _ can be found in [21| . 

Allowed magnetic dipole transitions between charmo- 
nium and bottomonium ground states have been consid- 
ered at NNLO in [22B2HI The results are: r(J/ip -> 7 r) c ) = 
(1.5 ± 1.0) keV and r(T(15) -> jr} b ) = (fc 7 /39 MeV) 3 
(2.50 ± 0.25) eV, where the errors account for uncertain- 
ties (which are large in the charmonium case) coming from 
higher-order corrections. The width _T( J /tp — > 7 rj c ) is con- 
sistent with 123]. Concerning 7^(T(1>S') — * 7%), a photon 
energy fc 7 = 39 MeV corresponds to a rjb mass of 9421 
MeV. 

The radiative transition Y(1S) 7I has been con- 
sidered in |25U26| . The agreement with the CLEO data of 
[2U is very satisfactory. J. Soto has reported about this at 
the conference. 



3 Low-lying QQq 

The SELEX collaboration at Fermilab reported evidence 
of five resonances that may possibly be identified with 
doubly charmed baryons [2E] ■ Although these findings have 
not been confirmed by other experiments (notably by FO- 
CUS, BELLE and BABAR) they have triggered a renewed 
theoretical interest in doubly heavy baryon systems. 

Low-lying QQq states are assumed to realize the hi- 
erarchy: m 3> mv 3> ^Iqcd, where mv is the typical in- 
verse distance between the two heavy quarks and vIqcd 
is the typical inverse distance between the centre-of-mass 
of the two heavy quarks and the light quark. At a scale \i 
such that mv ^Iqcd the effective degrees of free- 

dom are QQ states (in color antitriplet and sextet config- 
urations), low-energy gluons and light quarks. The most 
suitable EFT at that scale is a combination of pNRQCD 
and HQET |2^l3*D] . The hyperfine splittings of the dou- 
bly heavy baryon lowest states have been calculated at 
NLO in a s and at LO in ^Iqcd jm by relating them to the 
hyperfine splittings of the D and B mesons (this method 
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was first proposed in (HJJ). In [21], the obtained values are: 
- M- = 120 ± 40 MeV and Mc - Ma., = 34 ± 4 

"tc " cc "66 ~0» 

MeV, which are consistent with the quenched lattice deter- 
minations of 32,33,34,35 . Chiral corrections to the dou- 
bly heavy baryon masses, strong decay widths and elec- 
tromagnetic decay widths have been considered in [3*6], 

Also low-lying QQQ baryons can be studied in a weak 
coupling framework. Three quark states can combine in 
four color configurations: a singlet, two octets and a de- 
cuplet, which lead to a rather rich dynamics |2Sj. Masses 
of various QQQ ground states have been calculated with 
a variational method in |37j : since baryons made of three 
heavy quarks have not been discovered so far, it may be 
important for future searches to remark that the baryon 
masses turn our to be lower than those generally obtained 
in strong coupling analyses. 

4 High-lying QQ 

High-lying QQ states are assumed to realize the hier- 
archy: m 3> mn ~ Aqcb 3> rnv 2 . A first question is 
where the transition from low-lying to high-lying takes 
place. This is not obvious, because we cannot measure 
directly mv. Therefore, the answer can only be indirect 
and, so far, there is no clear agreement in the literature. 
A weak-coupling treatment for the lowest-lying bottomo- 
nium states (n = 1, n = 2 and also for the T(3S 1 )) appears 
to give positive results for the masses at NNLO in ^2] an d 
at N 3 LO* in 39 . The result is more ambiguous for the 
fine splittings of the bottomonium IP levels in the NLO 
analysis of and positive only for the Y(1S) state in the 
N 3 LO* analysis of 0Oj- In the weak-coupling regime, the 
magnctic-dipole hindered transition T(2S) — > 7 77b at lead- 
ing order |22] does not agree with the experimental upper 
bound while the ratios for different n of the radia- 
tive decay widths r(Y(nS) — > 7 X) are better consistent 
with the data if Y(1S) is assumed to be a weakly-coupled 
bound state and Y(2S) and Y(3S) strongly coupled ones 

Masses of high-lying quarkonia may be accessed by 
lattice calculations. A recent unquenched QCD determi- 
nation of the charmonium spectrum below the open flavor 
threshold with staggered sea quarks may be found in |43| . 
At present, bottomonium is too heavy to be implemented 
directly on the lattice. A solution is provided by NRQCD 
|4*4~] . Since the heavy-quark mass scale has been integrated 
out, for NRQCD on the lattice, it is sufficient to have a 
lattice spacing a as coarse as m 3> l/a> mv. A price 
to pay is that, by construction, the continuum limit can- 
not be reached. Another price to pay is that the NRQCD 
Lagrangian has to be supplemented by matching coeffi- 
cients calculated in lattice perturbation theory, which en- 
code the contributions from the heavy-mass energy modes 
that have been integrated out. A recent unquenched de- 
termination of the bottomonium spectrum with staggered 
sea quarks can be found in 03- Note that all match- 
ing coefficients of NRQCD on the lattice are taken at 
their tree-level value. This induces a systematic effect of 
order a s v 2 for the radial splittings and of order a s for 



the fine and hyperfinc splittings. In 03 > also the ratio 
r(Y(2S) - e+e-)/r(Y(lS) - e+e") x M* m /M* {1S) 
has been calculated. The result on the finest lattice com- 
pares well with the experimental one. 

In order to describe electromagnetic and hadronic in- 
clusive decay widths of heavy quarkonia, many NRQCD 
matrix elements are needed. The specific number depends 
on the order in v of the non-relativistic expansion to which 
the calculation is performed and on the power counting. 
At order mv 5 and within a conservative power counting, 
S- and P-wave electromagnetic and hadronic decay widths 
for bottomonia and charmonia below threshold depend on 
46 matrix elements |46| . More are needed at order mv 7 |47l 
48 , 49, . Order mv 7 corrections are particularly relevant for 
P-wave quarkonium decays, since they are numerically as 
large as NLO corrections in a s , which are known since long 
time |50| and to which the most recent data are sensitive 
|51llT] . NRQCD matrix elements may be fitted to the ex- 
perimental decay data [SlfllSH] or calculated on the lattice 
54,55,56 . The matrix elements of color-singlet operators 
can be related at leading order to the Schrodinger wave 
functions at the origin |3j and, hence, may be evaluated 
by means of potential models [S7| or potentials calculated 
on the lattice |5H|. However, most of the matrix elements 
remain poorly known or unknown. We refer to pQ for a 
summary of results. 

At a scale \i such that mv ~ Aqcd ^S> \i ^S> mv 2 , 
confinement sets in. Far from threshold, the effective de- 
grees of freedom are QQ states (in color singlet config- 
uration) and light quarks. Neglecting light quarks, the 
QQ propagation is simply described by a non-relativistic 
potential |59U6U| . This will be in general a complex val- 
ued function admixture of perturbative terms, inherited 
from NRQCD, which encode high-energy contributions, 
and non-perturbative ones. The latter may be expressed 
in terms of Wilson loops and, therefore, are well suited for 
lattice calculations. 

The real part of the potential has been one of the first 
quantities to be calculated on the lattice (for a review see 
[SHI)- In the last year, there has been some remarkable 
progress. In [ST], the 1/m potential has been calculated 
for the first time. The existence of this potential was first 
pointed out in the pNRQCD framework jSHj. A 1/m po- 
tential is typically missing in potential model calculations. 
The lattice result shows that the potential has al/r be- 
haviour, which, in the charmonium case, is of the same 
size as the 1/r Coulomb tail of the static potential and, 
in the bottomonium one, is about 25%. Therefore, if the 
1/m potential has to be considered part of the leading- 
order quarkonium potential together with the static one, 
as the pNRQCD power counting suggests and the lattice 
seems to show, then the leading-order quarkonium poten- 
tial would be, somewhat surprisingly, a flavor-dependent 
function. In |62| . spin-dependent potentials have been cal- 
culated with unprecedented precision. In the long range, 
they show, for the first time, deviations from the flux-tube 
picture of chromoelectric confinement |63| (for a review 
see, for instance, [OS])- The knowledge of the potentials in 
pNRQCD could provide an alternative to the direct deter- 
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mination of the spectrum in NRQCD lattice simulations: 
the quarkonium masses would be determined by solving 
the Schrodinger equation with the lattice potentials. The 
approach may present some advantages: the leading-order 
pNRQCD Lagrangian, differently from the NRQCD one, 
is renormalizable, the potentials are determined once for 
ever for all quarkonia, and the solution of the Schrodinger 
equation provides also the quarkonium wave functions, 
which enter in many quarkonium observables: decay 
widths, transitions, production cross-sections, ... . 

The imaginary part of the potential provides the NR- 
QCD decay matrix elements in pNRQCD. They typically 
factorize in a part, which is the wave function in the 
origin square (or its derivatives), and in a part which 
contains gluon tensor- field correlators [05,46 66 ,67 . This 
drastically reduces the number of non-perturbative pa- 
rameters needed; in pNRQCD, these are wave functions 
at the origin and universal gluon tensor-field correlators, 
which can be calculated on the lattice. Another approach 
may consist in determining the correlators on one set of 
data (e.g. in the charmonium sector) and use them to 
make predictions for another (e.g. in the bottomonium 
sector). Following this line in |65U68| . at NLO in ev s , but at 
leading order in the velocity expansion, it was predicted 
^had(X6o(2P))/r had ( Xb2 (2P)) m 4.0 and r had ( Xb i(2P))/ 
il la d(X62(2P)) ~ 0.50. Both determinations turned out to 
be consistent, within large errors, with the CLEO 111 data 
P- 

5 Threshold states 

For states near or above threshold a general systematic 
treatment does not exist so far. Also lattice calculations 
are inadequate. Most of the existing analyses rely on mod- 
els (e.g. the Cornell coupled channel model 6!3 or the 3 Po 
model E3). 

However, in some cases, one may develop an EFT ow- 
ing to special dynamical conditions. An example is the 
Jf (3872) discovered by BELLE |7T] and seen also by CDF 
|Z3, DO G3J and BABAR [Hj- If interpreted as a loosely 
bound D° D* and D° D* molecule, one may take advan- 
tage of the hierarchy of scales Aqcd ^> ^> m 2 / (2m re d) 
« 10 MeV > ^binding - Indeed, the binding energy, ^binding, 
which may be estimated from Mx"(3872) ~ (M D o * +M D o), is 
very close to zero, i.e. much smaller than the natural scale 
mj/(2m r0( j). Systems with a short-range interaction and a 
large scattering length have universal properties that may 
be exploited; in particular, production and decay ampli- 
tudes factorize in a short-range and a long-range part, 
where the latter only depends on one single parameter, 
the scattering length j?5Tl76ll77j . 

Another interesting case is provided by the F(4260), 
discovered by BABAR [23 and seen also by CLEO [23 
and BELLE [HOj- If interpreted as a heavy charmonium 
hybrid (see e.g. (EH)j one mav rely on the heavy-quark 
expansion and on lattice calculations to study its proper- 
ties. In particular, analogously to the energy of a static 
quark-antiquark pair, the energies of static hybrids have 
been calculated in quenched approximation on the lattice 



in [H2j- One may expect that the static energy is related 
to the static potential of the system and that this may be 
a relevant quantity for the dynamics of the system, but 
to substantiate this a suitable EFT for heavy hybrids, like 
pNRQCD for heavy quarkonium, needs to be formulated. 
Such a formulation does not exist yet. 1 



6 Production 

Although a formal proof of the NRQCD factorization for- 
mula for heavy quarkonium production has not yet been 
developed, NRQCD factorization has proved to be very 
successful to explain a large variety of quarkonium pro- 
duction processes (for a review see the production chap- 
ter in pP). In the last year, there has been a noteworthy 
progress toward an all order proof. In 83 , 84 , it has been 
shown that a necessary condition for factorization to hold 
at NNLO is that the conventional octet NRQCD produc- 
tion matrix elements must be redefined by incorporating 
Wilson lines that make them manifestly gauge invariant. 

Differently from decay processes, a pNRQCD treat- 
ment does not exist so far for quarkonium production. 
The difficulty of such a formulation may be linked to that 
of providing a full proof of factorization at the level of 
NRQCD and a consistent definition of the NRQCD pro- 
duction matrix elements at the level of pNRQCD. A pN- 
RQCD formulation of quarkonium production may present 
potentially the same advantages as that of quarkonium de- 
cay: a sensible reduction in the number of parameters and 
hence more predictive power. 

In the last years, two main problems have plagued 
our understanding of heavy quarkonium production: (1) 
double charmonium production in e + e _ collisions and (2) 
charmonium polarization at the Tevatron. 

In |85|. BELLE measures a(e + e~ — > J/ip+r] c ) Br(cc — > 
> 2 charged) = 25.6 ± 2.8 ± 3.4 fb and in [BE], BABAR 
finds er(e + e~ — > J/ip +r) c ) Br (cc — > > 2 charged) = 17.6 ± 
2.8^2'^ fb- When these cross sections first appeared, they 
were about one order of magnitude above theoretical ex- 
pectations. In the meantime, some errors have been cor- 
rected in some of the theoretical determinations, and, more 
important, NLO corrections in a s have been calculated in 
[8Tj and higher-order v 2 corrections in [HHj- All these im- 
provements have shifted the theoretical value much closer 
to the experimental one. In |89j . a preliminary estimate 
of a(e + e~ — > J/tp + r] c ) that includes the above correc- 
tions has been presented, it reads: 16.7 ± 4.2 fb. It ap- 
pears that, within the uncertainties, the discrepancy has 
been resolved. Still open is the issue of the inclusive double 
charm production in the presence of a J ~/ip. BELLE mea- 
sures a ratio er(e + e~ — » J/i/j + cc)/a(e + e~ — » J/tp + X) 

1 It is suggestive, however, that the lowest hybrid state is ex- 
pected to contain a pseudoscalar color-octet quark-antiquark 
pair and gluons, whose quantum numbers are those of the elec- 
tric cloud in a diatomic II U molecule, so that the system has 
J PC numbers like the Y. Moreover its mass, obtained 

by solving the Schrodinger equation, is consistent, within 100 
MeV, with the experimental range 4.25-4.30 GeV. 
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that is about 80%, to be compared with theoretical esti- 
mates, which are about 10% (see £Q for a detailed discus- 
sion). New theoretical analyses are timely. 

Charmonium polarization has been measured at the 
Tevatron by the CDF collaboration at run I on 110 pb _1 
of data [5U| and recently at run II on 800 pb _1 of data |§T] . 
The data of the two runs do not seem consistent with each 
other in the 7-12 GeV region of transverse momentum, px, 
and both are not with NRQCD expectations. For large px, 
NRQCD predicts that the main mechanism of charmo- 
nium production is via color-octet gluon fragmentation, 
the gluon is transversely polarized and most of the gluon 
polarization is expected to be transferred to the charmo- 
nium. For an analysis of the NRQCD prediction and its 
dependence on the adopted power counting, we refer to 
|92) . The CDF data do not show any sign of transverse 
polarization at large px- Before drawing definite conclu- 
sions, a polarization study from the DO experiment would 
be most welcome, at least to settle the possible discrep- 
ancy between the run I and run II data. 



7 Conclusions 

Many new data on heavy-quark bound states are provided 
in these years by the B-factories, CLEO, BES, HERA and 
the Tevatron experiments. Many more will come in the fu- 
ture from the LHC and GSI. They will show new (perhaps 
exotic) states, new production and decay mechanisms. 
What makes all this interesting is that we may investigate 
a wide range of heavy quarkonium observables in a con- 
trolled and systematic fashion and, therefore, learn about 
one of the most elusive sectors of the Standard Model: 
low-energy QCD. The tools for this systematic investiga- 
tion are provided by EFTs and lattice gauge theories. Still 
challenging remains for both the description of threshold 
states. 
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